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Abstract 

We consider the influence of non-equilibrium electronic neutrinos (and 
anti-ncutrinos) on the ncutron-to-proton ratio. These neutrinos would come 
from massive Vt annihilations Dt-Vt- — > f'gZ^e. For sufficiently large i't masses 
this new effect would strongly enhance the (n/p)-ratio, leading to a very 
stringent bound on the Vt- mass, even adopting a rather weak upper bound 
on the effective number on neutrino species during nucleosynthesis. 
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The tau-neutrino remains the only one which still can have mass in the MeV 
range. This possibility is theoretically viable, since it can be unstable on cos- 
mological time scales, thus avoiding the limits set by the relic density [|^. It is 
also quite interesting experimentally, due to the good prospects for improving the 
Vt mass limits at a B meson or tau-charm factory [^. In addition, an MeV tau 
neutrino can be quite interesting cosmologically Q. 

The present-day experimental limit on its mass is |^]: 

mu^<22,MeV (1) 
Primordial nucleosynthesis considerations rule out Vr masses in the range 

0.5 MeV < m^^ < 35 MeV (2) 

if ly-r is stable on the nucleosynthesis time scale and does not possess new inter- 
actions capable of depleting its density in the cosmic plasma 0. Calculations of 
spin-flip processes for the case of a Dirac-type Vr shift the lower bound in this 
interval down to 0.2 MeV [^]. Altogether these bounds can be summarised as 
m^^ < 0.2 MeV and mf,^ < 0.5 MeV, where "M" and "D" stands for the Dirac 
and Majorana mass respectively. 

All these bounds are obtained under assumption that the effective number of 
extra neutrino species Aki, during nucleosynthesis cannot be larger than 0.4 or 
0.6. Recent contradictory data on the primordial deuterium abundance ^] 
may cast some doubts on the validity of this limit (for recent analysis see refs. ||TT[| ). 
In particular, if Aki, = 1 is allowed, there may be an open window for neutrino 
mass somewhere near 20 MeV. 

In this letter we will show however that this is not the case, at least for a Dirac- 
type neutrino, and that high 1/^ mass values are excluded even when we adopt a 
much weaker upper bound on the extra number of neutrino species. Our result 
follows from the consideration of non-equilibrium v^s and Pe's originating from 
the annihilation of heavy i/^'s. The effect of non-equilibrium electronic neutrinos 
on the kinetics of (n-p)-reactions was considered in refs.|l^], for the case where 
the fe's and Pe's arise from heavy particle decays during nucleosynthesis. Non- 
equilibrium fg's and Pg's produced by e'''e~-annihilation in the standard model 



also change the (n/p)-ratio but the effect is small |13|. In contrast, the effect we 



^Of course such interactions are needed, at some level, in order to comply with the limit on 
the relic neutrino density, but the time scale involved is much larger. 
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consider here is much more significant. For example, for the case of a very heavy 
m ~ 20MeV Vr we show that it is equivalent to 0.8 extra massless neutrinos 
species for the case of a Dirac v-r-, and to 0.1 extra neutrinos for the Majorana 
case. The importance of this effect was mentioned in ref. but the calculations 
were not performed there because with the bound Afcj^ < 0.6 the complete range 
of high Vt mass values was already excluded even without this effect. 

One can understand the essential features of the relevant phenomena in the 
following way. At some stage v^s, being massive, went out of thermal equilibrium. 
Their number density tends to a "frozen" value, but their residual annihilation 
produced energetic non-equilibrium fe's and i>eS. If this take place sufficiently 
late, roughly speaking, when T < 2MeV, the produced non-equilibrium i/g and 
z^e are not thermalized. The resulting distortion of the z^e spectrum would shift 
the frozen (n/p)-ratio and, correspondingly, the primordial ^He abundance. The 
effect has different signs for high and low energy z^e's. If there is an excess of z^e's 
(or Pe's) at the high energy tail of the spectrum, the (n/p)-ratio increases. Indeed, 
this ratio is determined by the competition between the universe expansion rate 
and the rates of the reactions: 

p + e~ <-> n + z^e (3) 

n + e"*" ■«->■ p + (4) 

An excess of energetic z/g's and i^g's destroys neutrons in the first reaction and 
produces them in the second one, with the same rate. However, since the number 
density of protons at the moment of freezing of these reactions is about 6 times 
larger than that of neutrons, the second process is more efficient and we get more 
neutrons than in the standard case. On the other hand if there is an excess 
of neutrinos at low energy, its effect on the process (Q) is small because of the 
threshold effect. However, it enhances the destruction of neutrons in the reaction 
(Q) and, correspondingly, the (n/p)-ratio goes down. There is also another effect 
related to the increased overall number density of electronic neutrinos coming 
from the z^t-z^,- — > z^e^^e annihilation. This effect decreases the (n/p)-freezing tem- 
perature as well as the (n/p)-ratio independently of the energies of the electronic 
neutrinos. However, as we will see in what follows, the spectrum distortion at 
high energies is more important and the net result of all these effects is to enhance 
the (n/p)-ratio leading to more stringent z^,- mass limits. 

Though we obtained our results numerically, an accurate analytical treatment 
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is possible for large values of the Vt mass > lOMeV. The kinetic equation 
which governs the spectrum of electronic neutrinos has the form: 
dfei „ dfei 1 f (fpi d^P2 d^k 

Hid 



dt ^ duji 2wi J {2TT)nEi (27r)32S2 (27r)32a;2 

{2T:f 5\pi +p^-ki-k2)Y.\A^\ ifrlfr2 " /el/e2) (5) 

where and /g are the i^r and v^. spectral functions (or occupation numbers), 
respectively, and H is the Hubble parameter. If the expansion is dominated 
by photons, e~, e~^, and two massless neutrino species, the effective number of 
degrees of freedom is g^: = 9 so that H ~ 5T'^ /mpi. As a first approximation we 
take into account only the reaction i^ri + i^T2 ^ i^ei + i^e2 in the collision integral. 
Here A is the amplitude of this reaction and the summation is made over the 
spins of all particles except for v^i . For the case of a Dirac z^r one has 

Y^lAl \=32Glipik,f, (6) 

while for the Majorana case: 

I Al \= WGUiPih? + {P2kif - mHkik^)] (7) 

where the identical particle factor 1/2! has already been included in (0). 

The other relevant processes temporarily omitted in the right hand side (r.h.s.) 
of eq.(||), namely elastic scattering of Uf> (Pg) on themselves and other light 
fermions are essential at higher temperatures. They smooth down spectral dis- 
tortions and force the distribution back into equilibrium, somewhat diminishing 
the effect. We will take them into account in what follows (see eqs.(^, 10)). 

We will solve equation (||) (or equivalently eqs.(|^, |l^)) perturbatively assum- 
ing that /e has the form /g = flf^^ , where for the equilibrium part we take the 
Boltzmann expression flf^"^ = exp(— w/T). Quantum statistics corrections may 
diminish the effect by (5-10)% |lj]. We neglect terms of the second order in 5f . 
The distortion of the electronic neutrino spectrum is induced by the deviation of 
the massive Vr from thermal equilibrium which arises when the temperature of 
the cosmic plasma drops below the Vr mass. It is usually assumed that kinetic 
equilibrium is maintained, while an effective chemical potential of the same mag- 
nitude for particles and antiparticles is generated |15|. In other words takes 



the form = exp(^(t) — E/T). The Vr number density is calculated from the 
well known equation: 

rir + 3Hnr = {av) {n^^"^^ - nl) (8) 



where {av) is the thermaUy averaged cross-section of z^T--annihilation (multiphed 
by velocity) and nf'^'^ is their equihbrium number density. We write n-j- as usual, 
in the form n,- = rno where uq ~ Q.QQggT^ is the equilibrium number density of 
massless neutrinos with the same temperature and gg is the spin counting factor. 
We have calculated r numerically solving eq.(^) and parametrised the solution in 
the form r = tq + riuiu^/T. The results we obtain agree with those of ref. 

We assume that the temperature drops in the usual way, T = —HT, possible 
corrections to this law are not essential here. In this case the equilibrium part of 
/e annihilates the r.h.s. of eq.(^) and it would contain only 5f. Introducing new 
variables x = m/T and y = uj/T we can rewrite equation (^) in the form: 



(9) 

for the Dirac case while, for the Majorana case we obtain 

(10) 

Here a = x'^/y and the dots denote contributions from elastic scattering which 
we have so far omitted. 

In the limit of high mass, i.e. for x ^ 1 and a ^ 1, the integrals can be 
performed approximately so that the r.h.s. of these equations can be written as: 

{x - yf 



— ^ ^ 0.96 X lO-^m^ exp 

ox x^y^l^ 

(Dirac case), 

(x - yf 



FDix,y) + ... (11) 



— ^ ~ 1-92 X lO-^m^— ^ — 7A5-exp 
8x x°y~^''^ 



FMix,y) + ... (12) 



(Majorana case). 

Here and in what follows m is measured in MeV, Fd = (1 + lly/4a:^ + 
65y'^/32x^)/{l + 15/4x + 165/32x2), and Fm = (1 + 5y/ix^ - 35y^/32x^)/il + 
15/4x + 165/32x2). The numerators come from the expansion of the integrals in 



the r.h.s. of equations (^, 10) and the denominators come from the non-relativistic 
expansion of the number density. 

It is straightforward to take into account the other processes where and 
v'e participate, denoted by in the above equations. They are either elastic 

scattering off light leptons, fe + ^ ^ i^e + I, or i^e + z^e ^ I + I annihilations. 
These processes tend to restore equilibrium, thus diminishing the distortion 6f. 



The colhsion integral in the r.h.s. of eq. @ contains 5f in terms of two types: 
first, when 5f = 5f{ki) refers to the particle under consideration and, second, 
when 6f = 5f{k2) refers to one over whose momentum we integrate in d?k2. 
The first ones come with a negative sign and force the system into equilibrium, 
while the terms of the second type typically appear with positive sign and induce a 
deviation from equilibrium ^ In order to account for the terms of the second kind 
one has to solve a complicated integro-differential kinetic equation. In contrast 
it is very easy to treat the terms of the first kind. They result in the addition of 
an extra term of the form {—5f)F{x,y), with a known function F{x,y), to the 
r.h.s. of eqs.(|9,10). In what follows we take into account only the contribution of 
the terms of the first kind, thus somewhat overestimating the suppression factor. 
The real effect should be somewhat larger. However, even in this approximation 
our effect results in a considerable improvement of the strength of the bound on 

Taking into account all possible reactions of with light leptons we get 
F ~ m^y/3x^. The result is valid for arbitrary x and y values (the restriction 
X ^ 1 is not assumed here) . After a simple algebra we get the following equations 
for the spectrum distortion: 

where "M" and "D" stand for the Dirac and Majorana cases, respectively. 

The solutions of these equations can be easily written. These expressions 
should be substituted into the equation which governs the neutron-to-proton 
ratio. Since the spectrum is distorted at high neutrino energies we may neglect 
the electron mass, leading to the approximate form (see e.g. ref.[|l3[): 

dr„. 



dT 



(24 + 12^ + 2/32)(e-'3 - r„(l + e"^)) + 
(1/2) / dyy\y + Pf{5f{y + /3) - rn{5f{y) + 5f{y + /?)))! (15) 





Here is the ratio of the neutron number density to the total baryonic number 

density. The (n/p)-ratio is expressed through (n/p) = r„/ (1— r„), and (3 = Am/T 
*The terms of the second kind may also appear with a negative sign, like e.g. in the process 
1/^ + i/^ ^ i/e + i^e from neutrinos in the initial state, but the contribution from similar terms 
coming from neutrinos in the final state over-weights the negative one. 



5 



where Am = 1.3 MeV is the neutron-proton mass difference. Of course, in our 
numerical integration we have used the exact equation for the (n/p)-ratio as given 
e.g. in ref . . 

We can now roughly estimate the variation of (n/p)-ratio using eq.([l5|) in 
the following way. First, we integrate eqs. (13,14|) neglecting the last term in 



the r.h.s. In this approximation the t'e's and Pg's are in thermal equilibrium 
above T ~ 2 MeV, while below this temperature the elastic scattering processes 
are not important. We estimate the relevant integral over x in the Gaussian 
approximation near x = y and integrate the obtained result over y as written in 
eq.(l^). So this equation takes the form: 



^ = -O.OST^ 
dT 



(24 + 12/3 + 2(5^)e-^ + A-/2 - 
r„((24 + 12/3 + 2^2)(i + g-/?) + ^+/2 + A-/2)l (16) 



where 



A^ w O.OlTr^m^x ( 1 - — ± In — ) (17) 

^ X m XiJ 



Ai ^ OmArl,m' (in J ± 2/3 (1 - i)) (18) 

with Xi = m/Ti = m/2. The signs "+" and "— " in the above expressions 
correspond respectively to the integrals / dyy'^{y + P)'^Sf{y) and to / dyy'^iy + 
Pf6fiy + (3). 

As we have already noted, the modification of the neutrino spectrum results 
in two opposite effects. First, the "equilibrium" value of r„ (obtained by equating 
the r.h.s. to zero) increases. Second, the (n-p)-reactions get frozen later, as the 
excess of energetic i^e's and i^e's increases their rates. Taking into account both 
effects we can write the correction to r„ as: 



6r„, r. , A-eft-(A+ + A-)/(l + e-ft)T^^^_^^l + e^^_ 

1 + e-/3i ^ ^ 



1 



48 + 24/3i + 4/3f 

where Pq = Am/TfQ, Tjq = 0.65 MeV is the standard value of (n-p)-freezing 
temperature and A^ are taken at x = m/TfQ. The new freezing temperature Tfi 
can be found from 

Using these expressions we obtain 5rn/rn = (5 — 10)% for a Dirac Vr with mass in 
the range m = 10 — 20 MeV, while for the Majorana case we find 6rn/rn = (1 — 
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2)%. These results are in reasonable agreement with our numerical calculations 
described below. 



Now let us briefly describe the calculations. The solution of eqs.(13,14) was 
taken analytically in terms of one-dimensional integral Jq dx' . Then we substi- 
tuted this expression into the exact equation governing the (n/p)-ratio. Inte- 
gration over y in this equation can be approximately done analytically with a 
very good precision. The last integral over x' was performed numerically for 
different values of x and accurate interpolating functions were substituted into 
the equation for the (n/p)-ratio. This first order differential equation was solved 
numerically. 

In the case of Dirac Uj- we obtained Srn/rn = 4.2% for m = 10 MeV and 3.6% 
for m = 20 MeV. In the Majorana case we got 1.2% for 10 MeV and 0.5% for 
20 MeV. These results can be rewritten as corrections to the effective neutrino 
number Ngg as a function of the Uj- mass as shown in the figure. The solid curves 
denote the results we obtain including the effect of non-equilibrium electronic 
neutrinos and antineutrinos from massive Vr annihilations on the neutron-to- 
proton ratio. The dashed curves represent the corresponding results neglecting 
these corrections ^ • As can be seen from the figure the increase of N^g in the 
mass range of interest is larger for the case of Dirac than it is for the Majorana 
neutrino case. From our results one concludes that relatively long-lived massive 
tau neutrinos above the few MeV range are ruled out by the requirement N^g < 4. 
Thus it would seem from our results that the only way to accommodate a massive 
tau neutrino above the few MeV range is if there are new interactions beyond 
those of the standard model that can make it decay and/or annihilate efficiently 
on the nucleosynthesis time scale. The case of unstable z^r's has been considered 
in refs.[|l^]. The effect we have discussed would also improve the bounds obtained 
in these papers. The alternative case of neutrinos with large majoron annihilation 
cross sections is under consideration [ITi 



As this work was completed we became aware of the paper [18| where the non- 
equilibrium heating of electronic neutrinos by z^t-P,-- annihilation was considered. 
The conclusion of this paper about the influence of non-equilibrium fg's and z^g's 
on the (n/p)-ratio is opposite to ours. We think that the results differ because 
these authors took into account only a decrease of the (n-p)-freezing temperature 
due to the excess of electronic neutrinos, but neglected the spectrum distortion 
which we have discussed here and which turns out to have a stronger impact on 
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the (n/p)-ratio. 
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